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Abstract: Cubature formulas, asymptotically optimal with respect to accuracy, are derived for 
calculating multidimensional weakly singular integrals. They are used for developing a universal code 
for calculating capacitances of conductors of arbitrary shapes. 

1. Introduction 

Optimal with respect to accuracy methods for calculating singular integrals are being actively 
developed presently. They represent an important field of computational mathematics. Asymptotically 
optimal and optimal with respect to order ( to accuracy and to complexity) algorithms for calculating 
singular integrals on closed and open contours, and multidimensional singular integrals have been 
constructed in [1-3] on Holder and Sobolev classes of functions. 

In constructing optimal with respect to accuracy methods for calculating one-dimensional, bisin- 
gular and multidimensional singular integrals, a general method, proposed in monograph [1], was 
used. This method can be applied not only to singular integrals but also to other integrals with 
moving singularities. 
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This method allows one to construct several asymptotically optimal and optimal with respect to 
order and to accuracy algorithms for calculating hypersingular integrals [4], the Poisson and Cauchy 
type integrals [5], and multidimensional Cauchy type integrals [6]. 

Although multidimensional weakly singular integrals are used in many applications, optimal meth- 
ods for calculating these integrals are not developed. 

An exception is the book [1], where asymptotically optimal with respect to accuracy methods for 
calculating integrals of the form 
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< 71, 72 < 1, were constructed on Holder and Sobolev classes. 

Thus, the development of optimal methods for calculating multidimensional weakly singular inte- 
grals is an actual problem. Construction of efficient cubature formulas for calculating weakly singular 
integrals for calculating capacitances of conductors of arbitrary shapes by iterative methods proposed 
in [7] is very important in many applications, for example, in wave scattering by small bodies of ar- 
bitrary shapes and in antenna theory. A bibliography on methods for calculating capacitances and 
polarizability tensors is contained in [7]. 

In this paper the method proposed in [1] is generalized to multidimensional weakly singular 
integrals. As a result the analogs of the basic results for singular integrals, obtained earlier, are 
obtained for weakly singular integrals. Moreover, we study the applications of optimal with respect to 
order cubature formulas for calculating weakly singular integrals on Lyapunov surfaces. Our results 
are used for constructing an universal code for calculating capacitances and polarizability tensors of 
bodies of arbitrary shapes. 

This paper consists of two parts. 

In the first part of the paper optimal methods for calculating integrals of the types: 

Kf=JJ f(^^ida 2 < Sl)S2<27r . (L1) 
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are constructed on Holder and Sobolev classes of functions. 



Our results for integrals (1.1) can be generalized to the integrals with other periodic kernels and 
functions. The development of cubature formulas for integrals (1.1) is of considerable interest because 
the results are applicable to integrals with weakly singular kernels defined on closed Lyapunov surfaces. 
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It will be clear from our arguments, that the results can be generalized to I— dimensional integrals, 
Z = 3,4,-- - . 

The second part of this paper is deals with the iterative methods for calculating capacitances 
of conductors of arbitrary shapes. A general numerical method for calculating these capacitances is 
developed, and the results of numerical tests are given. 

2. Definitions of optimality. 

Various definitions of optimality of numerical methods and a detailed bibliography can be found 
in [8,9,10]. Let us recall the definitions of algorithms, optimal with respect to accuracy, for calculating 
weakly singular integrals. 

Consider the quadrature rule 

«i "-2 Pi P2 

T f=H J2J2J2P^ 2 hh(tut2)f {h ' h) (x kl ,y k2 )+ 
k 1= l k 2 =l h=0 h=0 

+R ni n 2 (/; Pk±k 2 hh ! x ki , Vk 2 '■> *1 > h), (2.1) 

where coefficients Pk lk2 hh (*i>*2) and nodes (x kl ,y k2 ) are arbitrary. Here f( ll,l2 \si, S2) = d ll+l2 f (s\, S2) / ds 1 ^ ds l 2 2 . 
The error of quadrature rule (2.1) is defined as 

Rn 1 n 2 (.f\Vk 1 k 2 hl 2 \ x k 1 ,yk 2 ) = SUp \Rn 1 n 2 (f;Pkik 2 l 1 l 2 ',X kl ,y k2 ;ti,t2)\. 

(tl,t 2 )6[-l,l] 2 

The error of quadrature rule (2.1) on the class ^ is defined as 

Rni^i^'iPhitehh'i x ki,yk 2 ) = SUp R ni n 2 {j \Pkik 2 hl 2 ',X kl ,y k2 ). 

Define a functional 

CnmaW = inf R nin2 (^;p kl k 2 hl 2 ;x kl ,y k2 ). 

Pk 1 k 2 l 1 l 2 ', x k 1 ,Uk 2 

The quadrature rule with the coefficients P klk2 i^ 2 an d the nodes (x^. ,y k ) is optimal, asymptot- 
ically optimal, optimal with respect to order on the class ^ among all quadrature rules of type (2.1) 
provided that: 

R ni n 2 i Pk^hh ' X *ki ' y*k 2 ) _ j lxl 

The symbol a x j3 means Aa < (3 < Ba, where < A, B < 00. 
Consider the quadrature rule 



Tf = 5>fc(ti,* 2 )/(M fc ) 



fc=i 

+R n (f;p k ;M k ;t 1 ,t 2 ), (2.2) 
where coefficients p k (ti,t 2 ) and nodes (-M&) are arbitrary. 



3 



The error of quadrature rule (2.2) is defined as 

R n (f;p k ;M k ) = sup \R n (f;p k ; M k ;t 1 ,t 2 )\. 

(tl,t 2 )6[-l,l] 2 

The error of quadrature rule (2.2) on the class ^ is defined as 

R n (V;p k ;M k ) = sup R n (f,p k ; M k ). 



Define a functional 



CnW = mf Rn(^; Pk ;M k ). 

Pk\M k 



The quadrature rule with the coefficients p* k and the nodes (M£) is optimal, asymptotically 
optimal, optimal with respect to order on the class \I> among all quadrature rules of type (2.2) provided 
that: 

R n (*;pl;M*) 

By i? nin2 (^ r ) the error of optimal cubature formulas on the class ^ is defined. It is obvious that 

Rn in2 {y) = Cmn 2 (*)- 

3. Classes of functions 

In this section, we list several classes of functions which are used below. Some definitions are 
from [11,12]. 

A function / is defined on A=[a,b] or A=K, where A' is a unit circle, satisfies the Holder conditions 
with constant M and exponent a, or belongs to the class H a (M),M > 0, < a < 1, if \f(x') — 
f(x")\ < M\x' - x"\ a for any x',x" G A. 

Class Hu, where u(h) is a modulus of continuity, consists of all functions / G C(A) with the 
property |/(xi) - f(x 2 )\ < Mu(\xi - x 2 \),x 1 ,x 2 G A. 

Class W r (M) consists of functions / G C(A) which have continuous derivatives /', /", . . . , /( r_1 ) 
on A, a piecewise derivative on A satisfying max x ^ a ^\f^'\x)\ < M. 

Class Wp(M), r = 1,2 . . . ,1 < p < oo, consists of functions f(t), defined on a segment [a,b] or 
one A = K, that have continuous derivatives /', /", . . . , f( r ~ 1 \ and an integrable derivative such 
that 



l/p 

< M. 



J \f (r) (x)\Pdx 
.A 

Class W£(M), r=l,2...,0<a<l, consists of functions /(t), defined on a segment [a,b] or 
one A = K, that have continuous derivatives /', /", . . . , f( r \ such that 

|/W(xi)-/W(x 2 )| <M|xi-x 2 | a . 

A function f(x±,x 2 , . . . , x/), Z = 2, 3, . . . , defined on A = [ai, 62; 02,^2! 
. . . ; ai, bi] or A = K\ x K 2 x • • • x Ki, where Aj, = 1, 2, . . . , I, are unit circles, satisfying Holder conditions 



with constant M and exponent ctj, i = 1, 2, . . . , I, or belongs to the class H aij ___ ;CXl (M), M > 0, < a < 
l,i = l,2,...,l, if 

|/(xi,x 2 , • • • - /G/1,1/2, • • • ,yi)\ < M(\ Xl - yi p +... + \ Xl - yi \ a <). 

Let u>, Ui, where i = 1, 2, . . . , I, I = 1, 2, . . . , be moduli of continuity. 

Class // Wl) ... ja , ! (M), consists of all functions / € C(A),^4 = [01,62502,62; • • • ;oz,&z] or A = K\ x 
K2 x ■ ■ ■ x Ki, with the property 

|/(xi,x 2 ,.. . ,x/) - f(y!,y 2 ,.. .,yi)\ < M{ui{\xi - yi\) -\ \-wi(\xi - yi\)). 

Let Hf(A),j = 1,2,3, A = [a 1 ,b 2 ;a 2 ,b 2 ; . . . ;a h bi] or A = K x x K 2 x ■ ■ ■ x K h I = 2, 3, . . . , be the 
class of functions f{x\,x 2 , . . . ,xi) defined on A and such that 

\f(x)-f(y)\ <co( Pj (x,y)),j = 1,2,3, 

where x = (xi, . . . , x/), y = (yi, ...,yi), pi (x, y) = maxi<i<i(\xi - yi\), p 2 (x,y) = Y!i=i N ~ 

yi\,p3(x,y) = ELi \ x i - yt\ 2 } 1/2 - 

Let Hf(A),j = 1,2, 3, A = [ai,b 2 ;a 2 ,b 2 ; . . . ; a h k] or A = K x x K 2 x ■ ■ ■ x K h I = 2, 3, . . . , be the 
class of functions f(x\,x 2 , . . . ,x{) defined on A and such that 

\f(x)-f(y)\<( Pj (x,y)r,j = 1,2,3. 

More general is the class Hpj(A),j = 1,2,3. It consists of all functions /(x) which can be rep- 
resented as /(x) = p(x)g(x), where g(x) € H?(A),j = 1,2,3, and p(x) is a nonnegative weight 
function. 

Let Z^{A),j = 1, 2, 3, be the class of functions f(x±, x 2 , . . . , x/) defined on A and satisfying 
\f(x) + f(y) - 2/((x + y)/2)\ < u{pj(x, y)/2),j = 1, 2, 3. 

Let Z"{A),j = 1, 2, 3, be the class of functions f(x\,x 2 , . . . ,xi) defined on A and satisfying 
|/(x) + f(y) - 2/((x + y)/2)| < ( Pj (x, y)/2) a ,3 = 1,2,3. 

Class Zpj(A),j = 1,2,3, consists of all functions /(x) which can be represented as /(x) = 
p{x)g{x), where g(x) G Z"(A),j = 1,2,3, and p(x) is a nonnegative weight function. 

Let W Tl '- ,ri (M), I = 1, 2, . . . , be the class of functions /(xi, X2, . . . , x/) defined on a domain A, 
which have continuous partial derivatives 

d\ v \f(x 1 ,...,x l )/dx^ ...dx?,0 < \v\ < r - l,\v\ = v ± + ■■■ + vi, ri > Vi > Q,i = 1,2, ... ,l,r = n + 
■ ■ - + ri and all piece-continuous derivatives of order r, satisfying \\d r f(xi, . . . , Xi)/dx r ^ . . . dx\ l \\c < M 
and \\d r <f(0, . . . , 0, Xi, 0, . . . , 0)/dx[ l || c < M, i = 1, . . . , I. 
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Let Wp 1 ''"' 1 " 1 (M), I = 1, 2, . . . , 1 < p < oo be the class of functions f(x±,x 2 , . . . ,xi), denned on a 
domain A = [ai, &i; . . . ; ai, b{\, with continuous partial derivatives 

<9H/( Xl) ... ^ti/dxl 1 ...dx^',0 < \v\ < r - 1, \v\ = vi + ■ ■ ■ + v h n > vi > 0,i = 1,2, . . . ,l,r = 
r± H + ri , and all derivatives of order r, satisfying 

\\d r f( Xl , Xl )/dxYdxZ . . . dx? \\ LM) < M, 

Wl+V2+ ... +Vlf{xiA Q)/dx r ldx v 2 dx v t \\ Lp{[aiM]) < M , \v 2 \ + H + ■ ■ ■ + h| < r - n - 1; 



ll^ 1+ - + ^- 1+r 7(o, . . . , o, Xl )/dx^dx% . . . dxl^dx? \\ Lp([aiM) < m, \ Vl \+\v 2 \+- ■ < r-n^-i. 

Let be A = [a±, 62; a 2 , b 2 ; ... ; a;, bi] or ^4 = K\ x x ■ ■ ■ x K\. Let C r {M) be the class of functions 
f(x\,x 2 , . . . ,xi) which are defined in A and which have continuous partial derivatives of order r. Partial 
derivatives of order r satisfy the conditions 

,,d\ v \f(x 1 ,...,x l ) l 



vi \\C 



< M 



11 dxi 1 ...Ox 1 ; 1 

for any v = (vi, . . . , w/), where > 0, i = 1, 2, . . . , I are integer and Yli=i v i = r - 
By ^ we denote the set of periodic functions of the class 

It is known [13] that Lyapunov spheres are defined as regions bounded by a finite number of 
closed surfaces satisfying the three Lyapunov conditions: 

1. At each point of the surface a tangent plane (and, therefore, a normal) exist. 

2. If 6 is the angle between the normals at the points mi and m 2 , and r is the distance between 
these points, then 

6 < Ar x , < A < 1, 

where A and A are positive numbers which do not depend on mi and m 2 . 

3. For all points of the surface, a number d > exists such that there is exactly one point at 
which a straight line, parallel to the normal at the surface point m, intersects the surface inside a 
sphere of radius d centered at m. 

Let S be a Lyapunov sphere, and N be the external normal to this sphere. We introduce a local 
system of Cartesian coordinates (XiVi 0) whose origin is located at an arbitrary point tuq of S, the ( 
axis is directed along the normal iVo at the point mo, and the x an d f? axes lie in the tangential plane. 
Within a sufficiently small neighborhood of mo, the equation of the surface S in the local coordinates 
(X) Vi C) nas the form 

( = F( X ,r))- 

Definition 4.1. [13] The surface S belongs to the class Lk(B,a) if F(x,r]) G W^(B), and the 
constants B and a do not depend on the choice of the point mo- 



4. Auxiliary statements. 

We need the following known facts from the theory of quadrature and cubature formulas. These 
facts can be found, for example, in [11], [12], [14], [15]. 

Lemma 4.1. Let *i be the class of functions W£(l), 1, 2, ... , 1 < p < oo, < t < 1, f(t) G 
and the quadrature rule 

\ n 

/ f(t)dt = Y,Pkf(tk) + Rn(f) 
fc=1 

is exact on all the polynomials of order up to p — 1, and has error Rn^i) on the class ^i. Let ^ be 
the class of functions W£(l), r = 1, 2, . . . , 1 < p < oo, a < x < b, and g(x) G W£(l). Then the 
quadrature rule 

b . n 

/ g(x)dx = ^p k g{a + (6 - a)t k ) + R n {g) 
a *=i 
/ias error R n {^2) on the class of functions ^2 and 

R n (*2) = (b-ay+ 1 - 1 /PR n (y 1 ). 



Theorem 4.1. [11] Among quadrature formulas 



/m p 
f( X )dx = ^ j>m/ (0 (^) + R w = L ^ + R (f) 



k=l 1=0 



the best formula for the class Wp(l) (1 < p < 00) with p = r — 1 and r = 1,2, • • • , or p = r — 2 and 
r = 2, 4, 6, • • • , is i/ie unique formula defined by the following nodes x* k and coefficients p* kl : 



:l = h(2(k-l) + [R rq (l)] 1 / r ), k = 1,2,... 



m, 



2h 2v+1 ( 
(Z = 0,l,...,p), p* kt2v = — ^"^(l), fc = 2,3,...,m-l; v = 0,l,... 



r — 1 



P*k,2v+i = I A; = 2,3, . . . ,m - 1; u = 0, 1, . . . , 



r - 2 



, h = 2 _1 (m - 1 + [R rq {l[ 



ll/rN-l 



r-1 



and R rq (t) is the Chebyshev polynomial t r + ^ deviating least from zero in the norm L q (—1, 1), 

i=0 

where p~ l + q -1 = 1. Here 



uw;(i)} = R n [w;{i) 



Rrq(l) 



2 r r\fFq~TT(m - 1 + [^(l)] 1 /' 



Let a function f(x, y) be given on a rectangle -D = [a, b; c, d]. Consider the cubature formula 

r. „ m n 

J J f(x,y)dxdy = ^2^2p ki f(x k ,yi) + R mn (f), (4.1) 



fe=i i=i 



defined by a vector (X, Y, P) of a nodes a < x\ < x 2 < • • • < x m < b, c < y\ < y 2 < • • • < y n < d, and 
coefficients pki- 

Theorem 4.2 [11]. Among all quadrature formulas of the form of (4-1) the formula 

,. ,. m n 

/ / f(x, y)dxdy = 4hqY,Yl f( a + ( 2fc " c + ( 2i ~ !)</) + ^(/), 
J J fe=l i=l 

where h = q = is optimal on the classes H^^^D) and H^{D). In addition 

h q 

Rmn[H^ 2 {D)\ =4mn[q J ^(fydt + h j uj 2 (t)dt}; 


q h 

R mn [H%(D)} =4mn J J uj^t 2 + T 2 )dtdr. 
o o 

Consider the cubature formulas of the form: 

JV 

J J p(x,y)f(x,y)dxdy = J> fc /(M fc ) + R(f), (4.2) 

D k=1 

where p(x,y) is a nonnegative and bounded on D function, p k , M k (M k G D) are coefficients and 
nodes. 

Theorem 4.3 [11]. Let p(x,y) be a nonnegative bounded weight function. If R]y[Hp j(D)] and 
R^[Zpj(D)], where j = 1,2,3, and < a < 1, are the errors of optimal formulas as (4-2) on the 
classes - {D) and Z'* j(D), respectively, then 



lim N a > 2 R N [H? AD)} = lim N a / 2 2R N [Z? AD)] = 

N^oo P ' J N^oo ' ' ' ■' 

- (2+a)/a 

j j(p(x,y)) 2 /^dxdy 



, j = 1,2,3, 



/ \ (2+a)/a 7r / 6 

where D l = ^ (^) / D 2 = 2 1 "«/(2 + «), and D 3 = 2 1 -/ 2 /(2 + a). 

7/j = 2, then the conclusion holds for n— dimensional cubature formulas. 
Remark. Theorem 4-2 is generalized to the case of unbounded weights p(x , y) in [2]. 
In this paper we will use the result of S.A. Smolyak (see [16]): 

Lemma 4.4 . Let H be a linear metric space, F be a convex centrally symmetric set 
with center of symmetry 9 at the origin, and L(f), l\ (/),..., In(J), be some linear function- 
als. Let S(li(f), . . . , Ijy (/)) be some method for calculating the functional L(f) using functionals 
(h(f), • • • , ^at(/)), and S be the set of all such methods. Then the numbers Di, . . . , D^ exist such that 

N 

sup \L(f) - V D k l k (f)\ = inf sup \L(f) - S(h(f), l N (f))\. (4.3) 
feF ^ s feF 



This means that among the best methods of calculating functional L(f) 

L(f)^S(h(f),...,l N (f)) (4.4) 

there is a linear method. 

Proof. Let us associate with each / G F a point (L(f), h(f), ■ ■ ■ ,In{/))- Let Y be a set of all 
such points (y , . . . , y N ) for / G F. 

From our assumptions, it follows that Y is a closed centrally symmetric set with the center of 
symmetry at the origin. 

Let 

D = sup y , 

( yo ,o,...,o)eY 

and consider the case Dq / and Dq / oo. 

Draw the support plane for the set Y through the point (Do, 0, . . . , 0) : 

N 

C (y -D ) + J2C j y j = 0. 
3=1 

One can choose a support plane for which Co / 0. Since Y is centrally symmetric with respect to the 
origin, the plane 

N 

C (y + D ) + J2c jyj = 

is also a support plane for Y, and Y lies between these two planes. 
Hence, we have for the points of Y the inequality: 

N 

\y - D W\ < D o, Dj = -Cj/C . 

3=1 

The definition of yi implies 

N 

sup|L(/)- V^(/)|<-Do- (4.5) 

Let /o be an element F corresponding the point (Do, 0, . . . , 0). Then S(h(±fo), ■ ■ ■ , Zjv(±/o)) = 
S(0, . . . , 0), and L(fo) — L(—fo) = 2Dq. Therefore either for fo or for — fo the error in (4.3) is not less 
than Do- This and (4.4) imply that the right-hand side in (4.3) is not less that the left-hand one. But 
the right-hand side of (4.3) can not be less than the left-hand side of (4.3) because a set of methods 
S for calculating functional (4.3) contains a set of linear methods. Lemma 4.4 is proved. ■ 

Corollary. Among all functions for which the optimal method for calculating L(t) has the greatest 
error for a given set of functionals, there exists a function satisfying the conditions l± (f) = • • • = 
^N(f) = 0. It follows from the proof that such a function is the function fo. 

5. Optimal methods for calculating integrals of the form (1.1). 



5.1. Lower bounds for the functionals Q nm and £jv. 

In this Section we derive lower bounds for the functionals Q n m and (n, denned in Section 2 , for 
calculating integrals (1.1) by the cubature formulas 

ni n 2 pi P2 

Kf= J2 Yl ^2^2pkik 2 hi 2 {si,s2)f ilul2 \x kl ,x k2 )+ 

k- L =l fc 2 =l «i=0 h=0 

+Rn 1 n 2 (f',Pk 1 k2hl2'i x ki,Xk 2 ;Si,S2), (5.1) 

and 

N 

Kf = Y.Pk(si,s 2 )f{M k ) + R N (f; Pk ; M k ; s u s 2 ) (5.2) 



k=i 

on Holder and Sobolev classes. 

Theorem 5.1. Let ^ = H UljW2 (D) or ^ = H^(D), and calculate integral (1.1) by formula (5.1) 
with pi = p 2 = 0. Then the inequality 

h q 



Cram 2 [*] > -^n x n 2 \q J ui(t)dt + h J u 2 (i)dt], 



where q = , h = , and 
^ ri2 ' "l ' 

2?r 2?r 

f f dsids2 , 

7 ' J J ( S m 2 ( Sl /2) + S m2( S2 /2))A 1 ' > 

o o 

is valid. 

Corollary. Let ^ = H aa (D) or ^ = H^(D), and calculate integral (1.1) by formula (5.1) with 
n\ = n 2 = n and pi = p 2 = 0. Then the inequality 

(1 + a)n a 

is valid. 

Proof of Theorem 5.1. Denote by tp(si,S2) a nonnegative function belonging to the class 
Hu lcJ2 (l) an d vanishing at the nodes (x kl ,x k2 ), 1 < k\ < n±, 1 < k 2 < n 2 - 
One has: 

Lin\U2 (V'i Vk\ki ! ^fci j 2- £12) — 

2tt 2tt / 2tt 2tt 

> J_ /" /" / /" /" ip(ai,a 2 )daida2 ■ , , 

- 4W y ly y a m2 ((<71 _ 5l)/2)+(5in 2 (((T2 _ 52)/2)) > i 

\o 

2tt 2tt / 2tt 2tt 

= 4^ / / ^1.^) / / -2^ _ .^^^(^ _ S2)/2)) A ) ^ = 

\o 

2tt 2tt 2tt 2tt 

1 /* f ds\ds 2 



4ir 2 J J (sin 2 (s 1 /2) + sin 2 (s 2 /2)) x 

00 00 



J J ip(s 1 ,s 2 )ds 1 ds 2 . (5.3) 
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From Lemma 4.4 and Theorem 4.2 one concludes that the following inequality 



T 

Rnm 2 (ip;Pkik2;x kl ,x k2 ) > -^n x n 2 



n q 

q J ui(i)dt + h J u 2 (t)dt 



IT IT 

,h = — , q = — 
ni n 2 



o o 
holds for arbitrary weights p kl k 2 an d nodes (x kl ,x k2 ) and 

h q 

Cnn(*) > ^n x n 2 \q j uj^dt + h J uj 2 (t)dt]. 



Theorem 5.1 is proved. ■ 

Theorem 5.2. Let ^ = Hf or Vl/ = Zf, i = 1,2,3, and calculate the integral Kf by cubature 
formula (5.2). Then 

( N [Hf] = 2( N [Zf] = (1 + (l)) 7 (47T 2 ) 2 / a AiV- a/2 , 

12 / 1 \(«+2)/2^ 6 d(p n _ 2 _ jn _ 2W2 



w/iere £>i = ^ J / ^fe^, D 2 = and D 3 = 

Proof. The proof of Theorem 5.2 is similar to the proof of Theorem 5.1, with some difference is 

2w 2tt 

in the estimation of the integral f f ip(si, s 2 )ds\ds 2 , where the function ip(si, s 2 ) belongs to the class 



Hf (or Zf), is nonnegative in the domain D = [0, 2ir] 2 , and vanishes at N nodes M k , k = 1, 2, . . . , N. 
Using Lemma 4.4 and Theorem 4.3, one checks that the inequalities 

2tt 2tt 

inf sup f [ i/j(s 1 ,s 2 )ds 1 ds 2 = (1 + o(l)) A(4^ 2 ) (2+a)/o iV~ Q/2 , 

M k xpeHf ,ip(M k )=0 J J 
2w 2n 

inf sup f j i){si,s 2 )ds l ds 2 = {l + o{l))]-D i {AiT 2 ) i2+ol V a N- a l 2 

M k tbez? MM k )=o J J 2 
00 

hold for arbitrary M k € D, k = 1, 2, . . . , N. 

Substituting these values into inequality (5.3), we complete the proof of Theorem. ■ 
Theorem 5.3. Let ^ = C' 2 (l), and calculate the integral Kf by formula (5.1) with p\ = p 2 = 0, 

and n\ = n 2 = n. Then 

U[*]>(l + o(l))^, 

where K r is the Favard constant. 
Proof. Let 

1p(si,S 2 ) = V'l(si) + Y>2(S2), 

where < ipi(s) G W r (l) vanishes at the nodes x k , k = 1, 2, . . . , n, and < ip 2 (s) G VF r (l) vanishes 
at the nodes y k , k = 1,2, ... ,n. 

According to [11], for arbitrary nodes x k , k = 1, 2, . . . , n one has: 

2w 

J Ms)ds>—^, i = l,2. 
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Thus, the inequality 

2tt 2tt 

/ / ip(si,S2)dsids2 > —~ r — 



holds for arbitrary nodes (x±, . . . , x n ) and (j/i, . . . , y n ). 

The conclusion of Theorem 5.3 follows from this inequality and from (5.3). ■ 

Theorem 5.4. Let ^ = Wp T (l), r = 1, 2, . . . , 1 < p < oo, and calculate the integral Kf by 

formula (5.1) with p\ = p2 = r — 1 and n\ = ri2 = n. Then 

where R rq {t) is a polynomial of degree r, least deviating from zero in L q {[— 1, 1]). 

Proof. Let L = [j^]. Take an additional set of nodes £ fc = ^r, fe, / = 0, 1, . . . , L — 1. 

By (vi,Wj), i,j = 0, 1, . . . , iV — 1, iV = n + L, denote the union of the sets (xk,yi) and (£«,£?)• 
Let i/j(si,S2) = Y>i(si) + V'2(s2), where V ; i( s ) £ ^p(l) vanishes with its derivatives up to the order 
r — 1 at the nodes Uj, i = 0, 1, . . . , N — 1, and tp2{s) G Wp^(l) vanishes with its derivatives up to 

Vi + l 

order r — 1 at the nodes Wj, j = 0, 1, . . . , N — 1. Assume J i/ji(s)ds > 0, i = 0, 1, . . . , N — 1, and 

Wj+1 

J %l)2{s)ds > 0, j = 0, 1, . . . , N — 1, where vn = 2ir and wn = 27r. 



Let 

0, if (<7i, cr 2 ) = (si,s 2 ), 

-, otherwise, 



h(si,S 2 ,(7l,(72) := 



^ (si,s 2 ) = < 
^~(si,s 2 ) = 



(sm 2 ((ui-si)/2)W((ff2-S2)/2)) A 

^(Sl,^), ^(«1,S2)>0, 

0, i/ ^(si,s 2 )<0, 



0, i/ ^(si,s 2 )>0, 
-t/j(s 1 ,s 2 ), if tp(s 1 ,s 2 ) < 0. 
For each value i, j = 0, 1, . . . , N — 1, we have: 

2tt 2tt 

h(^,i,^j , o-i, o- 2 )ip((Ti, a 2 )daida 2 = 



If 





AT-1 jV-1 &+ 1 



^2^2 J J h (Ci^j^ cr i^ cr 2)ip(o-i,o- 2 )do 1 da2 



k=0 1=0 t c 



N-1N-1 5fc + l5 '+ 1 

^2^2 / ^ ^ cr 2)V' + (o-i, o- 2 )do 1 da 2 - 



fc=0 «=0 , c 



N-1N-1 ffc + l5i + 1 

^2^2 / h (^i^j^ cr i^ cr 2)ip~(o-i,o-2)da 1 da2 > 



fc=0 Z=0 t t 
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i+[(JVi-l)/2] j+[(JV 2 -l)/2] 



> ^2 h (Ci,Cj,Ck+i,Cl+i) / / 4> + (cr 1 ,a2)da 1 da2 + 

k= l+ i l= j+ i I I 

i+[(JVi-l)/2] j-\ J 

+ fr(&>€?>£fc+l,&-l) / / V' + (cri,0-2)rfo-i ( icr2 + 

fc=i+l i=j-[(tf 2 -l)/2] / fc ^ 

i-1 j+[(JV 2 -l)/2] €* ft+i 

+ X] X] M&>£.?>£fc-i>£H-i) / / ^ + (cr ll a 2 )da 1 da 2 + 

fc=i-[(JVi-l)/2] i=j+l £ 

i-1 i-1 ft? ft 

+ X] M&>€?>£fc-l>£/-l) / / ^ + (cri,cr 2 )(icricZc72- 

fc=i-[(JVi-l)/2] i=j-[(JV„-l)/2] J_ i 

i+[(JVi-l)/2] j+[(JV 2 -l)/2] ffc+i ft+l 



k=i+l l=j+l 



i+[(JVi-l)/2] j-1 ^+1 ft 

^2 M&>£i>£fc>6) / / ip~{cri,(T2)daid(72- 

fc=i+l J=j-[(JV a -l)/2] / fc ^ 

i-1 i+[(AT 2 -l)/2] ft? ft+i 

M&>€?>£fc>6) / / V 7 " (o-i, o- 2 )do-i dcr 2 - 

fc=i-[(JVi-l)/2] l=j+l £ 



•i-i i-i 



X] ^(&> 6) J J ip (a 1 ,a 2 )da 1 da2 



fc=i-[(JVi-l)/2] ;=j-[(JV 2 -l)/2] 



6-1 Si-i 



i+[(jV!-i)/2] j+[(jv 2 -i)/2] ftn-i ft+i 

= X] / / ip(cri,cr 2 )da 1 da 2 + 

k= l+ i l=j+i I I 

i+[(JVi-l)/2] i-1 ft=+l ft 

+ X] X] fr(&>€?>£fc+l>&-l) / / Y^l, cr 2 )(icri(io- 2 + 

fc=i+l l=j-[(JV a -l)/2] ( k ^ 

i-1 i+[(iV 2 -l)/2] ft? ft+i 

+ X] X] M&>£i>6-l,6+l) / / ^(«7i,cr 2 )dc7idcr 2 + 

fe=i-[(iVi-l)/2] I=j+l I 

i-i i-i % I 

+ ^(Ci, Cj, Cfc-i, 6-1) / / ^(cri,(72)dcrido-2- 

fc=i-[(iVi-l)/2] J=j_[(JV a -l)/2] f /_ x ^ 

i+[(JVi-l)/2] j+[(JV 2 -l)/2] 

(M6>^&>£0-fc(6>^&+i>6+i))>< 

fc=i+i /=i+i 

?fe+i 6+ 

x 



y y ip (cr 1 ,a 2 )da 1 da 2 - 
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i+m-i)/2] j-i 

12 (M&> - M&>&>&+1>6-1))> 
k=i+l l=j-[(N 2 -l)/2] 

£i 

£k+i £1-1 

i-l j+[(JV 2 -l)/2] 

X] (M&>&>&>6)-M6>&>&-i>6+i)): 

fe=i-[(ATi-l)/2] Z=j+1 



y y ^ (cr 1 ,a 2 )da 1 da 2 - 



i-i 



h(&, 6-1)): 

fc=i-[(JVi-l)/2] i=j-[(AT 2 -l)/2] 



x y y V' (fi, (J 2 )do\do 2 = 

£k-l £l-l 

= Ji + h + -h + Ja + h + h + h + h- 



Let us estimate the integral 



£k £l 



£k + l £l + l 



J J ip {a\,a 2 )da\da 2 < j J \ip (ai,a 2 )\daidcr 2 < 



£k £l 



< 



£k+l £l + l £k+l £l + l 

y y i^(o-i,o- 2 )idcrid<72 < (6+1-6) y 1^)1^ +(6+1-60 y iv^i^ < 



6 £; 



/27r\ r+2 1 
- 2 UJ H' 

where we have used the fact that the functions ipi(s) and ip 2 (s) on the segments 6+1] an d [6)6+1] 
vanish with derivatives up to order r — 1. 
Now let us estimate the sum: 

i+[(L-l)/2] j+[(L-l)/2] 

X] X] IM6,£j,6fc+i,6+i) - M6>£?;6fc,6)l = 
fc=i+i 1=3+1 



i+[(L-l)/2]j+[(L-l)/2] 

E E 

fc=i+l i=j+l 



sin z — 



(k+l-i) 



K(l+l-j) 



< 
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< 



i+[(L-l)/2] j+[(L-l)/2] 

E E 



;m /J 

i+[(L-l)/2]j+[(L-l)/2] 

C 
I 



fc =i+l i= i+ i ( sm 2!H +Sffi 2!(H)j 

^2+2A 



1+A 



/c — z Z — 7 
+ 



L L 



< 



(k-i) + (l-j) 



L 



< 



/[(£-l)/2] n [(i-l)/2] n . 

<« a E m+ E i £ 



<<:(£) 



2A 



=1 fc=l 

' ^ 2A i/ A < I 
logL i/ A = |, 
( 1 i/ A>±. 
By c > various estimation constants are denoted. Thus 

1 



L =0 



rv 



The expressions 12,13, and I4 are esimated similarly. 

From the definition of the function V'( s i) s 2) it follows that the error of cubature formula (5.1) 
for si = £i, S2 = can be estimated as follows: 

2w 2tt 

R(i>,,ti,£j) = J J Hvi,V2)h(£i,£ j ,a 1 ,a 2 )da 1 da2>o(J^J + 





i+[(L-l)/2]j+[(L-l)/2] 



k^+i J=i+ i / fc I 

i+[(L-l)/2] j-1 & 

+ X] X &+l>6-l) / / ^(<7l,<7 2 )<^1^2 + 

i-1 j+[(L-l)/2] €k &+i 

+ E E h (€h€j,€k-i,€l+i) / / ^(cri,a 2 )da 1 da 2 + 

k=i-[(L-l)/2] l= 3+ l ll 

i-i j-i Hi 

+ E E h (€h€j,€k-l,€l-l) / / 1p(<Tl,<T2)d<7ld(T2- 

k=i-[(L-l)/2]l=j-[(L-l)/2] "' ' 



Averaging the above inequality over i and j, one gets: 



Rnn[®} > SUp max -RnnO, > 



j L-1L-1 

^ 71 EE 



j=0 j=0 



i+[(L-l)/2]j+[(L-l)/2] 

E 

k=i+l l=j+l 



E E ^(&>&>6+1>£h-i) y y Y>(0"l,O"2)dc7icfcr2+ 
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i+[(L-l)/2] j-1 ^+1 & 

+ X] X ^(&>£.?>£fc+l>6-l) / / 1p(<Tl,(T2)daid(72 + 

k= l+ l i=i-[(L-l)/2] / fe ^ 

i-1 j+[(L-l)/2] & &+i 

fc=i-[(L-l)/2] J=i+1 ^ £ 

i— 1 j — 1 5i 

+ X X h (Ci,Cj,Ck-l,Cl-l) / / l/)((Tl,<T2)d(7ld<T2 +0 [~f) 

fc=i-[(L-l)/2] i=j-[(L-l)/2] {fc-ift-i 



> 



1 

L2 



L-l L-l 



i+[(L-l)/2]j+[(L-l)/2] ffc+i6+i 

„■ i 1 / ,■ ii ^ ^ „■ — n „■ — n 



5fc ft 



i+[(L-l)/2] J=j-1 ^+1 «i 



i=o i=o 



L-l L-l 



+ X / / ^( <J l' (7 2)^ 1 d(72^^/l(^,Cj,6+l,6-l) + 



fc=*+l j-[(L-l)/2] ?fc ^ 



i=0 j=0 



L-l L-l 



i-1 j+[(L-l)/2] & _ 

X X / / ^1,^2)^1^2 ^^M&,£j,£fc-l,6+l) + 



fc=i-[(L-l)/2] i=i+l {fc _ 1 fi 



i-1 



ft ft 



i=0 j=0 
L-l L-l 



+ ^2 ^2 i } {^i,^2)da 1 da 2 ^2^2h(Ci,Cj^k-i,Ci-i) = 

k=i-[(L-l)/2] i= J -[(L-l)/2] 4 - / _ i ^ i=0 i=0 



2vr 2?r 



°(~2~) + ^~2 y y 1p(c r l,C r 2)d(Tld(J2X 







2?r 2?r 



do\do2 



v0 



(smVi/2) + sm 2 (a 2 /2)) ; 



+ 



logn \ 
n ) 



where the following relation was used: 

L-l L-l 



^£!>&,^a-i,6-i) = o( 

i— n o— n V 



logn 



n 



2tt 2tt 



+ 



da±da 2 



[sm 2 ((Ji/2) +sm 2 ((T 2 /2)] 



A ' 



(5.4) 



i=o j=o 

Without loss of generality one may assume k = 1, Z = 1 in the previous equation. Let us estimate 

2tt 2tt 



t/o 



2 L-l L-l 



L 2 



i=0 j=0 



do\d(J2 







(sin 2 (a 1 /2) + sin 2 (a 2 /2)) 



< 



< 



L-l L-l - . i ; 
i=0 ^° ft I 



(sin 2 mm+sin 2 ((Q/2)Y 
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(sin 2 (a 1 /2) + sin 2 (a 2 /2)Y 



do\da 2 



+ 



6 6 



o o 



(sin 2 (a 1 /2) + sin 2 {a 2 /2)) 



-do\do 2 



+ 



ui + u 2 , 



where Y^' means summation over / (0,0). 
Let us estimate u\ and u 2 . One has: 



L-1L-1 

EE'/ / 



i=Q j=0 



(sm 2 ((ai)/2) + S m 2 ((a 2 )/2))' 



+ 



{sin 2 {i i /2) + sin 2 {a 2 /2)) x 

L-1L-1 



do\do 2 



+ 



EE' 

i=0 j=0 



1 



K((( ! )/2) + S m 2 (( ( r 2 )/2))' 



(«n 2 (6/2) + «n 2 ((,-/2))' 



da d a 2 



The expressions uu and ui 2 can be estimated similarly. Let us estimate u\\: 

L L 



nil 



<-,yy 

L toU K(te)/2) + ,m 2 (fe)/2)) 



T+A — 



L L 



< 



EE' 



< C- 



T2-2X (j2 i ? -2\l+A — t2-2\* 

i=0 i=0 V J ' 



where c > stands for various estimation constants. Hence 



U\ < 



L 2-2\ ■ 



Let us estimate u 2 : 



u 2 







(sin 2 {a 1 /2) + sin 2 (a 2 /2)) 



-da d G 2 



< 



< c 



Using polar coordinates, one gets: 



// (a 2 + , 
o o v 1 



-do\da 2 . 
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Thus: 



U < 



L 2-2\- 

From Lemmas 4.4, 4.1, and Theorem 4.1 it follows that 

2tt 



/ 



i)i{ai)dai > 



(l + (l))(2vrr+ 1 /^ rg (l) 
2 r r!(rg + l) 1 /?^ - 1 + [R rq (l)]^y ' 



where R rq (t) is a polynomial of degree r, least deviating from zero in L q ([— 1, 1]). 
Theorem 5.4 follows from inequalities (5.4) and (5.5). 



(5.5) 



5.2. Optimal cubature formulas for calculating integrals as (1.1). 
Holder class of functions. 

Let x k := 2kir/n, k = 0, 1, . . . , n, A w = [x k , x k +i,xi,xi +1 ], k, I = 0, 1, . . . , n - 1, x' k = (x k+ i + 
x k )/2, k = 0, 1, . . . ,n - 1, and (si,s 2 ) G Ay, i, j = 0, 1, . . . , n - 1. 
Calculate the integral if / by the formula: 

*/ = EEV*,*0 / / : daid ° 2 : — + *">■ ( 5 - 6 ) 



fe=0 «=0 



2 y 1 a "'" \ 2 

Theorem 5.5. Let \& = H aa (D),0 < a < 1. Among all cubature formulas (5.1) with pi = P2 = 0, 
formula (5.6), which has the error 

(2 + o(l))7 /7T\« 



i?T7 



(2 + OWj7 /7TX« 

ll J i + a w ' 



is asymptotically optimal. Here 7 is defined in (5.1'). 

Proof. Using the periodicity of the integrand, we estimate the error of cubature formula (5.6) 
as follows: 

/ (0-1,0-2) - /(ar-,x'- 



|-Rnn| — 



n—1 n—1 

EE 

fc=0 1=0 



si n 



2£r-£i 



2£2--£2 



< 



<( 

EE 

fc=0 1=0 



szn^ — 2~ l - + Sl1l z 2" 



+ s«n 



da\da2 



< 



n—1 n—1 



f (0-1,0-2) - f{x' k ,x\) 
A fc! v 2 2 ' 



-do\do^ 



+ 



+ 



n—1 n—1 



EE / / (/(«)- /(44)) : 

fe =o z=o 4 ; 



sin 



) ^sin 



+ sin 2 



,0-2-X, 



do\doi 
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= r\ + r 2 . 

Let us estimate each of the sums r\ and r2 separately. One has: 



n < 



+ 



i+M j+M 

E E 

k=i-M l=j-M A 



(siu 2 ^ + sin 2 



-da±da 2 



+ 



n—l n—l 

EE' 

fc=0 z=o 



( sin 2£l=2i +8in 2Za=S2.y 



-do\do 2 



r\\ + ri2, 



where EJ EJ' means summation over (k,l) such that 
A M ^ A*, A* = [xi- M ,Xi+M+i;xj-M,Xj+M+i],M = [Inn]. 
Furthermore 

da\du 2 



m < — 



< 



A* ^ 



sin 



2wM/n 2tt 



< 







dpdtp clogn 

p2\—l j^a-\-2—2X y 



Estimating ri2, one can assume without loss of generality = (0,0), and get: 

ir/nir/n 



n—l n—l 



/ (uJi(ai) +uj 2 {^2))daida 2 '^2'^2h M (si,S2,cri,a2) < 





n/nn/n 



k=0 1=0 



n—l n—l 







71- n 2+a 



/n II, ± II, ± 

/ {of + )daida 2 • s 2 ) o 1: a 2 ) < 

J k=0 1=0 

n—l n—l 

^2^2h kl (s 1 ,s 2 ,a 1 ,a 2 ) < 

+ 0{L) /TT\ a I j d(Jld<J 2 

1 + a \n) 



1 + a \n> 

1 + 0(1) /TP 



Here 



hki(si,s 2 ;ai,a 2 ) = sup s 2 ; ai, a 2 ). 

(o-i,o- 2 )eA fci 

Combining the estimates of r n and r 12 , one gets: 

l + o(l) /tt\« 

H < ; ; — 2 ( - ) 7 



1 + a Vn> 

Let us estimate r 2 . To this end we estimate the difference 



r 2 (M) = / / |/(4,^)"/(^,4)|x 
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sin z — k + str- 



da±da2 



First, we estimate 



»"2(*,j) < 



(an 2 (T1 2 Sl + sin 



2 o~2~S2 N ^ 
2 



» ( 



■sm z — 2— l - + sin A — 2 



< 



2+o-2A ' 



n 



The value r2(/c, X) is estimated similarly for \k — i\ < 3 and \l — j\ < 3. 
Let us estimate r2(&, Z) for other values of k and I. 
One has: 

r 2 (fc,0 = / / |/(4,^)- 704^)1 x 



> ( 



szn z — 2— 1 + szn z — 



do\doi < 



\k~i\\ + f\l-j\ 



n 



n 



1 



2 ai-x' i +d 1 (s 1 -x' i ) 



+ sin 2 ^^ 



l+A 



sin 



< 



c 

< — 

n 6 



\k-i\ 



n 



+ 



n 



\k-i\ 



n 



+ 



\1~J\ 



n 



n 



\k - i\ 2 + \l - j\ 2 



l+A 



< 



< c (| fc _ i | + |f_ J -|)i+ a 



- n o+2-2A /,, -i 2 



;|2^1+A - 



(\k-i\ 2 + \l-j\ 2 ) 
_c (|fc-z| 2 + |/-j| 2 )( 1+Q )/ 2 < 



— n a+2-2A 



C 



— n a+2-2A 



{\k-l\ 2 + \l-j\ 2 ) 
1 



(|jfe + j^-a/T+X ■ 



To estimate r 2 , one sums up the last expression over k and /. Without loss of generality assume 
= (0,0). Then 

/ [n/2]+l [n/2]+l 

^^2T 16 + 4 E E ^ 



1 



n L 



(k 2 + /2)A+l/2-a/2 I ' 



fc=0 Z=0 

where ^ J2' means summation over k and / such that k > 3 or / > 3. 
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One has: 



[n/2]+l [n/2]+l 

E E ' 

fc=0 «=o 



(&2 + /2)A+l/2-a/2 



< 



< A 



[n/l]+l [n/2]+l [n/2]+l 

X/ Z»2A+l-a ~*~ E E 



fc=3 



fc=3 «=3 



(fe2 + ;2)A+l/2-a/2 



< 



< A i 



1, i/ 2A - q > 1; 
logn, i/ 2A — a = 1; 



1-2A-R* 



i/ 2A-a<l. 



Hence 



n -(a+2-2A) ) if 2A-a>l; 
r 2 <-4^ n-Mogn, i/ 2A - a = 1; 

n -1 , i/ 2A - a < 1. 



Thus, if a < 1, then 



r 2 < o(n^). 
Combining the estimates of n and r 2 , one gets: 

«.Hls 7- (2 + 0(1)) 



1 + q \n> 

Theorem 5.5 follows from the comparison of this inequality with the lower bound of the value 
Cnn[H a ,a{D)], mentioned in the Corollary to Theorem 5.1. ■ 

Remark. If a = 1, the cubature formula (5.6) is optimal with respect to order. 

The proof of Theorem 5.5 yields also the following result: 

Theorem 5.5'. Let = H aa (D),0 < a < 1. Among all possible cubature formulas (5.1) with 
Pi = P2 = 0, formula 

*/ = EE/044) / / : — I 1 2 l +flnn, 



fe=0 «=0 



(sm 2 (^)+sm 2 (^))' 



which has the error 



Rnn[V] 



(2 + o(l))7 /7T 



1 + a \n 
is asymptotically optimal. 

To apply formula (5.6), one has to calculate the integrals 

da\du2 



i-kl 



IS 



O CT1 — X. 
l2_i L 



) CT2-X 



(5.7) 



A fcl ^-^ + ^ 2 

for A;, Z = 0, 1, . . . , n — 1. Exact values of these integrals for arbitrary values A are apparently unknown. 
Therefore the procedure of numerical calculation of integrals (5.7) should be given for practical appli- 
cation of formula (5.6). 
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Let k = i and / = j. Then the integral Jy is replaced by the integral 

7T 7T 

n n 

f f da 1 da 2 
Pij = / / t , h > 0, 

(sm 2 ^+sin 2 f ) A + /j 

which can be calculated by cubature formulas ( in particular, Gauss quadrature rule ) with arbitrary 

degree of accuracy because the function 
1 

(sin 2 ^+sin 2 ^-) 

has derivatives up to arbitrary order. The choice of parameter h is discussed in Section 8. 
Let k = i,l j, and 

Let k i, I = j, and 

47T 2 / . 2 x' k -x'A~ X „ 

/fc J = "^-( sm 2^) =p *i- 



Let k j£ i, I j£ j , and 



4vr 2 / . 9 x', - x- . ,xi - a;' 



-A 



The integral Kf is calculated by the formula 

n—1 n—1 

K ^ = EE P ^^'^ + Rnn{f,Pll,X k ,y'i). (5.8) 

Formula (5.8) is not optimal since it is not exact on constant functions f(x,y) = const. But one 
can estimate the error of this formula: 

n—1 n—1 

where M = max \f(x, y)\. 

The values — p* kl \ are easily estimated, and one gets the conclusion of Theorem 5.5'. 
Classes of smooth functions 

Theorem 5.6. Assume (p £ W r,r (l). Let ^ = W r,r (l), and calculate the integral Kip by formula 
(5.1) with p\ = r — 1, p2 = r — 1, and n\ = n 2 = n. Then the cubature formula 

2tt 2tt 

J J (stn 2 (a l - si)/2 + sm l [o 2 - s 2 )/2) A 



is asymptotically optimal. 

Before proving Theorem 5.6, let us describe the construction of the spline ip m n- Let xu = 
2/c7r/n, k = 0, 1, . . . , n. Divide the sides of the squares f2 = [0, 2ir; 0, 2ir] into n equal parts. Denote by 
A*,; the rectangle Am = [2kir/n, 2(fc+l)7r/n; 2l7r/n, 2(Z+l)7r/n], fc, I = 0, 1, . . . , n — 1. Let (si, s 2 ) E A y -. 
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First we approximate <p(<Ji, 02) as a function of 02, and construct a spline </? n (0i, 2 ) by the following 
rule. Let 01 be an arbitrary fixed number, < 01 < 27r. On the segments [2k7r/n, 2{k + 1)tt/ti] for 
fc ^ j — 2, . . . , j + 1, one has: 



Vn(01,02) = ^ 



Z=0 



y(°'')(0i,2fc7r/n) 



(02 - 2k7r/n) 1 + BiS^(a u (k + 1)/ 



7) 



where 



xr \ f \ ^(p^ l \a u 2kir/n) t 
(5(01, 02) := 99(01, 02) - 2^ j\ (°"2 - 2/cvr/n) 

1=0 



The coefficients are defined by the equation 

r— 1 

(2(* + 1W« " - E r-T-TTJ- (M* + l)/n - 02)^^ = 

(r — / — lj! n 

= (-l) r i? r i (2vr(2fc + l)/2n; vr/n; 2 ) , 

where R r i(a,h,x) is a polynomial of degree r, least deviating from zero in the norm of the space L 
on the segment [a — h,a + h}. On the segment [27r(j — 2)/n, 27r(j + 2)/n] the function ^(01,02) is 
defined by the partial sum of the Taylor series: 

, c\ • 1 \ y {Q > l) {<Ji,2iTj/n) . . . , 
¥M0i,02j = ^(0i,2vrj/n) H - (02 - j/n) H h 

H ; tv; (o"2 - 27rj/n) 



(r-1)! 

We define the function </? nn (0i, 02) by analogy with the function <p n ((Ti, 02)- 

Proof of Theorem 5.6. Let (51,^2) £ A^-. The error of formula (5.9) we estimate by the 
inequality 

n-l n-1 , , 

Vl "!) C2j - <Ain(01, 02 J 



\Rnn\ ^yZ^Z 



k=0 1=0 



(- 



sin 2 + sin 2 



-da\&02 



+ 



-1 n— 1 



^(01, 02) - 9 ? nn(01, 02) 



n + r 2 , 



(5.10) 



fc =0 /=o (sin 2 2if^ + sin 2 2*=2i) 

Ei 
means summation over (&, Z) such that i — 1 ^ fe ^ « + 1, ^ Z ^ n — 1 or < ^ n - 1, 



En 
means summation over the other values of (k, I). 



k,l 



Let us estimate each of the sums n and r 2 separately. In addition without loss of generality 
assume that ff (99(01, 02) — <Pnn(0i, 02))^0id02 ^ 0. Then 



i-l n-l 



n ^ XI XI l^^l' °"2) - Vnn(01, 02)1 J J 



do\do<i 



k=0 1=0 



A fe! ( 



,2 o-i-si 



-.2 cr 2 -A-2 



sin" " 1 2 ,J1 + sir 2 
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< A < 



i+ l+[( n _l)/2] j+l + [( n _ l)/ 2 ] 

^4 E E 



n-( r+1 ) ,A<l/2 

n -(r+2-2A) )A> l/ 2 ; 



(5.11) 



fc=i+2 



(sh 



2 Xfc-si , - 2 a;-s 2 



i+l+[(n-l)/2] j+l+[( n -l)/2] 

- 4 E E 

fc=j+2 i=j+2 



=i+2 ( sin" + sm 

1 



J J ip(<Ti,a 2 )da 1 d(72- 



( s in 2 + sin 2 



) (sin 2 ^±pi + sin 2 x -^^ 



J J " ip (a 1 ,a 2 )da 1 da 2 = r 2 i+r 22 , 



(5.12) 



where ip(ai,a 2 ) = (p(ai,a 2 ) (^l,^) 

■0 + (Cl,^) = < 
^~ (0l,<72) = < 



■0(o-i,<t 2 ) , if i>(<Ji,a 2 ) > 

, if V(oi,o- 2 ) < 0; 

, ifV(<7i,<7 2 ) >0; 

-^(0-1,0-2) , if -0(0-1, o- 2 ) < 0. 



It is obvious 

j+l+[(n-l)/2] j+l+[(n-l)/2] 

1 E E , 

fc=i+2 /=j+2 (sin" + sin ^ 

Let us estimate the integral 



1 + 0(1) 
A - 4?r 2 



2tt 2tt 



da\da 2 







2 £1 _l_ cir, 2 £2 



(5.13) 



sin -4- + sm 



i = JJ ip(ai, a 2 )do\ da 2 



99(0-1,0-2) - ^ n (o-l,0- 2 )^do-icZo-2 



+ 



+ 



J j (^(01,0-2) ~ *Pnn(vi,CT2f)d(T 1 d(J-2 



h + i 2 - 



(5.14) 



Since the expressions i\ and i 2 are estimated similarly, we estimate only i\. One has: 



2tt 

i\ ^ — max 

n «i 



y ^(si,<T 2 ) - <Pn(si, 0-2)^0-2 



This integral is a continuous function of s±, which attains its maximum at a point s*, and 



2vr 



n 



Xl + l 



^ / M°'*V,"S)| 

r!n J 1 1 



(x i+ i - o- 2 ) 
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r-1 



Bij(x l+1 - x{)r\ 



i=o 



(r-l-i)! 



(zj+l - <7 2 ) 



r-j-1 



r!n 7 f-' (r - 1 - jj! 



— / |-Rrl(<7 2 )|(i(T2 



4 / tt \ r+2 



From inequalities (5.14) and (5.15) one gets: 



(r + 1)! \n 



^FTWln) ^(D 



(5.15) 



(r + 1) 



and 



7*21 ^ 



2tt 2tt 



(r + 



v sin 2 ^+sin 2 f 



A ' 



One has: 



r 22 = o(n r ). 
Estimate (5.17) follows from the inequalities: 



(5.16) 



(5.17) 



J J ip (a 1 ,a 2 )da 1 da 2 ^ f [ |#7i,<7 2 )|d<7idff2 = 0(n- r " 2 ) 



and 



i+l+[(n-l)/2] j+l+[(ra-l)/2] 

E E 

fc=i+2 i=j+2 



gin 2 x^-s, + gin 2 



) ( 



sin — - + sin - ±L 



1-^2 y 

2 J 



< An^ E ' :/ '' /j ' " J \ < < < 



The estimate 



r ((*-i) 2 +(i-j) 2 ) 



i?nnW < (1+0(1))- 



A+l 



,A< 1/2 
nlogn , A = 1/2 
,A> 1/2. 



n 



2A 



27T r i? r l(l) 



:7 



(r + l)!(n-l + [i? rl (l)]i/r)r 
follows from inequalities (5.10), (5.11), (5.16), and (5.17). 

Theorem 5.5 follows from the comparison of the values Cnn[^] and i? nn [^]. ■ 
Let us construct cubature formulas for calculating integrals Kf on classes of functions W rr (l). 
These formulas will be less accurate than the ones in Theorem 5.3, but they will be optimal with 
respect to order, and easier to apply. 
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First, we investigate the smooth function 

2tt 2tt 



f(ri,T 2 )dTidT 2 



J J ( sin 2n_!i + sin 2^) A 

v 2 2 ' 

assuming f(t±,t 2 ) S W r,r . Changing the variables T\ = T\ — t, t 2 = t 2 — t, in the last integral, one 
gets: 

2 TT 2 TT 

fin +h,r 2 + t 2 )dTidr 2 



*P(ti,t 2 ) = j f 







(sin 2 2- + sin 2 f) A 
Thus, i/j(h,t 2 ) e W r > r . 

Remark. It is known [9] that Kolmogorov and Babenko widths on the class of functions W r,r (l) 
are equal to 5 n (W r ' r (l)) x d n (W r ' r (l), C) x ^75- Hence the recovery of the function ip(ti,t 2 ) using n 
functional is not possible with accuracy greater than O (^—^72^ ■ More precise conclusions are obtained 
in Theorems 5.3 and 5.4- 

Thus, for recovery of a function ip(ti,t 2 ), (ti,t 2 ) £ [0,27r] 2 with the accuracy <9(n~ r / 2 ), it is 
sufficient to calculate the value of the function tp(ti,t 2 ) at the nodes (vk,vi), where = 2kir/N, 
k, I = 0, 1, . . . , N, and iV 2 = n, and to use the local spline ipN(ti,t 2 ) of degree r with respect to each 
variable. 

Let us describe the construction of such spline. 

Assume for simplicity that M := N/r is an integer, and cover the domain [0, 27r] 2 with the squares 
A-ki = [wk, wi], k,l = 0, 1, . . . , M — 1, here = 2kir/M, k = 0, . . . , M. Approximate the function 
ip(ti,t 2 ) in each domain by the interpolation polynomial ipiy(ti,t 2 , ^-h) constructed on the nodes 
(xf, x l j), i,j = 0,1,... ,r, x\ = w k + jfci, i = 0,1, ... ,r. 

Denote the local spline, which is defined by the polynomials ipN(ti, t 2 , Am), by ipN(ti,t 2 )- 
If the values tp(vk,V[) are calculated by formula (5.9) with the accuracy 0(ra" r / 2 ), then 

Il^(tl,t2)-^(tl,t2)||c7<0(^ r/2 )- 

Therefore the spline ipN(ti,t 2 ) is optimal with respect to order, and a method for recovery of the 
function ip(t\,t 2 ), which has the error 0(n~ r l 2 ) (in the sup —norm) is constructed. 

6. Optimal methods for calculating integrals of the form Tf. 

Lower bounds for the functionals ( mn and (jy. 

First we get a lower bound for the error of formula (2.1) with p\ = p 2 = and n\ = n 2 = n, on 
Holder classes. 

Theorem 6.1. Let ^ = H aa (D), and calculate the integral Tf by formula (2.1) with n\ = n 2 = n 
and pi = p 2 = 0. Then the estimate: 

1 1 



c m > (1 + 0(1)) I [ dhdt2 (en 

W J " (l + a)W 7 (r 2 + t 2 )A [bA) 
-1 -1 
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holds. 

Proof. Let n be an integer number, L = [n/logn]. Let v k '■= — 1 + 2k/ L, k = 0,1,..., L. 
By ((,k,rji) we denote a set which is the union of nodes (xi,yj), i,j = 1,2, ...,n of formula (2.1) 
and the nodes (vi,Vj), i,j = 1,2,..., L. Let A M = [v k ,v k +i;vi,vi +1 ], k,l = 0, 1,...,L- 1. Let 
< VK^i)^) G H aa (D), where D = [—1, l] 2 , vanishing at the nodes fe, Z = 0, 1, . . . , iV, where 

N = n + L. 

Consider the integral 

l l 

/ rp i \ / ^ f f ip(T 1 ,T 2 )dT 1 dv i 



((n-vi) 2 + (r 2 -^-) 2 ) A 
1 -1 



L-1L-1 L-lj-1 i-1 L-l j-1 i-r 

EE + EE+EE+EE 



// 



1p(Tl,T 2 )dTldT 2 

({ri-ViY + ^-v^Y ~ 
A fc; 



A fe+i,i+j 



+ EE(x) ((i+1)I + (l + 1)J)i y y «r„r 2 )dr 1(i r 2+ 



i-1 L-j-1 / . 2A 



fc=0 «=0 



+ E E ( 2) (tkTWTiTTTWJ J ^ Tl ' T2)dTldT2+ 



i-i i-i / L \ 2A 1 f f 

+ EE(2j (( ,+ 1)2 + (;+1) w y «»•-»)*>*>= 

^i-fc-lj-i-l 



^^fL\ 2X U(L-i-l-k)U(j-l-l) f f w J 

+ EEU WrizTTTTW / / ^,r 2 )^r 2+ 



(( fc+ l)2 + ( Z + 1 )2)A 

L-1L-1 , , N 2A 



fc=0 i=0 x 7 A 

^i-fc-l,j'+! 



+ EE(§) -[-^y / / — - 

A;_fc_ij_;_i 

Here 17(A;) = 1 for k > 0, and J7(fc) = for k < 0. 
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Averaging the above inequality over all i and j, i, j = 0, 1, . . . , L — 1, one gets: 

j L-l L-l 

Pki;x k ,yi) > -2 E ^T(^)(Ci,??j) > 



i=0 j=0 



L-l L-l 



> 



\2\X 



L 2-2A 2 2A EE,,...,, 

k=o i=o u 7 v ; 

xU(L-j-l-l)J J 4>(T 1 ,T 2 )dT 1 dT 2 + 



L-l L-l 
i=0 j=0 



L-l L-l 



+ E E ^ " * " 1 " k ) U( j ~ 1 " l ) J j i>{T U T 2 )dT l dT 2 + 



i=0 j=0 
L-l L-l 



ifc+i,j-i-i 



+ J2J2 U ( i - 1 - k M L -j- 1 - 1 ) j j ^(Tl,T 2 )dT 1 dT 2 + 



i=0 j=0 



L-l L-l 



H-k-l,j+l 



+ E E U <j - 1 - - 1 - / / ^(n, r 2 )dTidr 2 



i=o j=o 



L-l L-l 



&i-k-l,j-l-l 
1 1 



> 



-2Ao2A EE 



i 2 - 2A 2 2A fr- ^((fc + i) 2 + (i + i) 2 ) : 

fL-fc-2 1 

+ y J 1p(n,T 2 )dTldT 2 + J J 1p(Tl,T 2 )dTldT 2 + J J ^{j\,T 2 )dT\dr 2 



1 



^(tI) T2)dri(ir2+ 

VL-k-2 Vl-I-2 



Vk -1 



-1 «i 
L-l L-l 



> 



2Ao2A EE 



£2-2A 2 2A 



fc=0 2=0 

From inequality (6.2) it follows that 



((fc + l)2 + ( Z+1 )2)A 



L-l L-l 



1 -1 
1 1 

1p(T 1 ,T 2 )dTldT 2 . 

-1 -1 



1 1 



Cr, 



[H aa (D)} > (l + (i))_ ^_ E E (P + / 2) a / / ^(n,7a)d7id75 = 



l + o(l) / } _dhdh_ 

4 7 7 (t?+ti) A j 

-1-1 -1-1 



1p(.T 1 ,T 2 )dT 1 dT 2 . 



From Theorem 4.2 and Lemma 4.4 it follows that the inequality 

l l 



J J 1p(Tl,T 2 )dTldT 2 



> 



4 1 

1 + an a 



-i -i 



valid for an arbitrary vector of the weights and the nodes (X, Y, P) on the class H aa (D). 
Theorem 6.1 follows from inequalities (6.3) and (6.4). ■ 
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Theorem 6.2. Let $ = C' 2 (l), and calculate the integral Tf by formula (2.1) with p\ = P2 = 0. 
If n\ = ri2 = n, then 

l l 

ds\ds2 



-i -l 

where K r is the Favard constant. 
Proof. Let 

ip(si, s 2 ) = tpi(si) + ^2(^2), 

where < ipi(s) £ W r (l), vanishes at the nodes x^, k = 1, 2, . . . , n, and < ip2(s) G VF r (l) vanishes 
at the nodes yk, k = 1, 2, . . . , n. 

For arbitrary nodes x k , k = 1, 2, . . . , n, one has (see [11]): 

1 

/2K r 
ipi(s)ds > - — i = l,2. 
(7rn) r 

-1 

Thus the inequality 

1 1 

ip(si,s 2 )dsids2 > 



II 



r 



(7rn) r 



-1 -1 

holds for arbitrary nodes (x±, . . . , x n ) and (yi, . . . ,y n ). 

Theorem 6.2 follows from this estimate and inequality (6.3). ■ 

Theorem 6.3. Let ^ = Wp' r (l), r = 1, 2, . . . , 1 < p < 00, and calculate the integral Tf by 

formula (2.1) with p\ = p2 = r — 1 and n\ = n2 = n. Then the estimate 

1 1 

, r Tl , , « 2 1 /9 J R ro (l) f f dsids 2 . N 

* (1 + ° { » K,r q+ !)./.(„- 1+ PMIF^ J J WTW' <6 ' 5) 

-1 -1 

holds, where R rq (t) is a polynomial of degree r, least deviating from zero in L q {[—1, 1]). 

Proof. Let L = [n/logn]. Consider the nodes (vk,vi), Vk = ^7, k,l = 0,1,..., L — 1. By 
(CiiVj)j hJ = 0,1, ■ ■ ■ ,N — 1, N = n + L denote the union of the nodes (xk,yi) and Let 
i>(si, s 2 ) = i>i(si) + ip 2 (s 2 ), where < ipi(s) € Wp"(l) vanishes with its derivatives up to order r — 1 
at the nodes i = 0, 1, . . . , N — 1, and < ip2(s) G W£ (1) vanishes with its derivatives up to order 

Vi+l 

r — 1 at the nodes rjj, j = 0, 1, . . . , N — 1. Assume that f ip\(s)ds > 0, z = 0, 1, . . . , N — 1, and 

Vi 

j Ms)ds>0, j = 0,l,...,N-l. 

Wj 

Using the argument similar to the one in the proof of Theorem 6.1, one gets: 

^ L-1L-1 

Cnn(^,Pkf,V k ,Vi) > ^2^2T(lp)(vi,Vj) > 

i=0 j=0 

* Z^x g g WTwhiTWY I / ^ T2)dTldT2 = 



-1 -1 
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11 11 
1 + o(l) f f dtidt 2 



(tj + t^A 



y y 1p(Tl,T2)dT 1 dT 2 . (6.6) 



"-2 

-1-1 -1-1 

Prom Theorem 4.1 and lemma 4.4 it follows that the inequality 
l l 

JJ^rM > (1 + „(!)) - - ^ (6.7) 

-1 -1 

is valid for arbitrary weights and the nodes (X, Y, P) on the class H aa (D). 
Theorem 6.3 follows from inequalities (6.6)- (6.7). ■ 
Cubature formulas. 

Let us construct a cubature formula for calculating the integral Tf on the Holder class 
H aa (D). Let Xk := —1 + 2k/n, k = 0,1, ... ,n, x' k = (x^+i + x k )/2, k = 0, 1, . . . ,n — 1, and 
Am = [x k ,x k +i;xi,xi +1 ], k,l = 0, 1, . . . ,n - 1. 

Calculate the integral T/ by the formula 

n— 1 n— 1 

r/ = EE/04 *D / / ((ri _ tl) 2 T + ( ^ 2 _^A +^(/)- 



(Yti — ti 

fc=0 z=o ^ / 



Consider another cubature formula for calculating the integral Tf. 

Let (<i, £2) S Ajj. By A* denote the union of the square Ajj and of those squares Am which have 
common points with the Ay. Consider the formula 

dT\dri 



Tf = fix'^x'j, 

A* 

n—l n—1 



fe=0 «=0 



where ' means summation over the squares which do not belong to A* . 

Theorem 6.4. Among all cubature formulas (2.1) with p\ = p 2 = and n\ = n 2 = n, formula 
(6.8), with the error estimate (6.1), is asymptotically optimal. 

Remark. Similar statement holds for formula (6.9). 

Proof of the Theorem 6.4. Let us estimate errors of formulas (6.8) and (6.9). 
The error of formula (6.8) can be estimated as follows: 

|Pnn(/)|<EE / / //,_^U^-^2U ^ + 



((ri-ti) 2 + (r 2 -t 2 ) 2 ) A 
|/(r 1 ,r 2 )-/(x' fc ,x' / )| 



+ EE / / T7Z—T^r-TZ— -T^dndr, = n + r 2 . (6.10) 



k=0 1=0 A kl 



k=0 1=0 
n— In— 1 „ „ 

((ri-ti) 2 + (r 2 -t 2 ) 2 y 

where Y^' means summation over k and / such that the squares Am belong to A*, and YJ" means 
summation over the other squares. 
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Let us estimate r\ and r 2 : 
2 



n < 



dT\dr 2 



< 



A* 



((r 1 -t 1 ) 2 + (r 2 -t 2 ) 2 ) A " n 2 - 2A +« 



o(n- Q ); 



(6.11) 



n—l n— 1 



fc=o ;=o 

Here h(Aki) denotes the maximum value of the function ((n — t\) 2 + (r 2 — i 2 ) 2 )~ A in the square A^j. 
One has: 

1 



L((n-ti) 2 + (r 2 -t 2 ) 2 ) A 



h(A 



dT\d,T2 



II 



< 



[ ((ri - ti) 2 + (r 2 - t 2 ) 2 ) A ((x k - h) 2 + ( Xl - t 2 ) 2 ) x \ 
2X(x k -ti + qi(n - x fc ))(n - Xk)) 



II 



+ 



({xk -h + gi (n - x fc )) 2 + (r 2 - t 2 ) 2 )^+i 
2\(xi -t 2 + g 2 (r 2 - x/))(r 2 - x/)) 



dT\dr 2 



dridr 2 



+ 



< 



({x k - h) + (xi -t 2 + q 2 {r 2 - xi)) 2 ) x 



dT\dr 2 



< 



2A(ti - x fc ) 

((x k -h + q(n - Xk)) 2 + (r 2 - t 2 ) 2 )A+i/ 2 



dT\dT 2 



2A(r 2 - xQ) 

((ri - ti) 2 + (z, - t 2 + g 2 (r 2 - xO) 2 ) A+1 / 2 



dT\dr 2 



+ 



< 



2 4 A n 2A+1 



< - 



2 4 A 



( fe 2 +/ 2)A+l/2 n 2-2A ( fc 2 + /2)A+l/2 ' 



where it was assumed that k > i + 1, and / > j + 1. Estimates for the other combinations of k and I 
are similar. Thus: 

n—l n— 1 



dr\dr 2 



r2 "(l + «)« a ^"i 7 ((TL-flF + fa-ta) 2 )* 



+ 



2 3 A 



n—l n—l 

EE"; 



(l + a)n 2 - 2A + Q (p + /2)A+i/2- 

Let us estimate the last term in the previous inequality. 
One has: 

1 



(6.13) 



n/2 n/2 

t t ^ + ~ fc ^/ 2 ] ^5//^ +/2 ) A+1/2 



n—l rt— 1 

EE"; 



< 
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1, A>l/2 

<cl logn, A = 1/2 (6.14) 
^ n l ~ 2X , A < 1/2 
where Yl XT means summation over k and I, (k, I) / (0, 0). 

In deriving (6.14) we have used the known result ([17], Theorem 56) which says that a number 
of points with integer-value coordinates, situated in the circle x 2 + y 2 = r 2 , is equal to nr 2 + 0(r). 
From inequalities (6.13) and (6.14) it follows that 

(i + (i)) l r r d Tl dT 2 



k=0 1=0 

This and (6.11) yield: 



(1 + a) n«^£- 7 7 ((n-^ + ^-^A" 



l l 



Rnn[Haa{D) = — — SUp / / — ^ < 



l l 

1 + o(l) f f dr\dr2 



~ (l + a)n"J J (r 2 + r 2 y (fU5) 

-l -l 

Theorem 6.4 follows from a comparison the estimates of C[H a ^ a (D)] and R nn [H a ^ a (D)]. ■ 
Let us construct asymptotically optimal cubature formula for calculating integrals Tf on the 
classes W rr . In the derivation of formula (5.9) the local spline if n {h,t2), approximating the function 
(p(t±,t2) in the domain [0, 2ir; 0, 2tt], was constructed. A spline f n n(tijt2), approximating the function 
f(t\,t2) in the domain [—1, 1] x [—1, 1], can be constructed analogously. Calculate the integral Tf by 
the formula 



l l 

fnn{n,T 2 )dT 1 dT2 
((n - tl) 2 + (r 2 - t 2 f 



Tf= I I JJZ " ^2 ' /" + *-(/)■ ( 6 - 16 ) 



-1 -1 

Theorem 6.5. Let ^ = W r ' r (l),r = 1, 2, . . . , and calculate the integral Tf by formula (2.1) with 

pi = pi = r — 1, and m = ri2 = n. T/ien cubature formula (6.16), which has the error 

l i 



fl „„ w ^ (1+< , (1)) n ^';> /r)r // 



( r + l)!(„-l + [fl rl (l)]lA)^ 7 (Tf + rJ)*' 

-1 -1 

is asymptotically optimal. Here R rq (t) is a polynomial of degree r, least deviating from zero in 
L q ([-l,l}). 

As in the proof of the Theorem 5.6 one gets the following estimate 

l l 

2R rl (l) f f dndT 2 



-i-i 

Comparing this estimate with the estimate of ( nn [W r,r (l)], from Theorem 6.3, one finishes the 
proof. ■ 
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7. Calculation of weakly singular integrals on piecewise continuous surfaces. 

In Sections 5 and 6 of the paper asymptotically optimal methods for calculating weakly singular 
integrals defined on the squares [0, 2ir] 2 or [—1, l] 2 were constructed. 

It is of interest to study optimal methods for calculating weakly singular integrals on piecewise- 
Lyapunov surfaces. 

Consider the integral 

Jf= [f- /(Tl ' T2 ' T3)dS r T ,ti,t 2 ,t3€G > (7.1) 



G ((Ti-ti) 2 + (r 2 -t 2 ) 2 + (T3-t 3 ) 2 )' 



where G is a Lyapunov surface of class L s (B,a). 

We show that the results derived in Sections 5 and 6 can be partially generalized to the integrals 
(7.1). 

Calculate integrals (7.1) by the formula: 

n p 

J f = E E Pkvf (v) (M k ) + R n (f,G,M k ,p kv ,t), (7.2) 

k=l \v\=0 

where t = {h,t 2 ,t 3 ), v = (v 1: V2,v 3 ), \v\ = vi + v 2 + v 3 , f ( - v \h,t 2 ,t 3 ) = q^t^q^ ■ 
The error of formula (7.2) is: 

R n (f,G,M k ,p kv ) = sup\R n (f,G,M k ,p kv ,t)\. 

tec 

Assume / £ Vl/i, and G € ^f 2 . Then the error of formula (7.2) on the classes and ^2 is: 
R n (t> 1 ,V 2 ,M k ,P kv ) = sup R n (f,G,M k , Pkv ). 



Let 



Cn[*i,*2] := s inf R n {^l,^2,M k , Pkv ). 

M k ,p kv 



A cubature formula with nodes M£ and weights p* kv is called optimal, asymptotically optimal, 
optimal with respect to order on the class of functions *f/i and surfaces ^2) if 

R n ^ l ^ 2 ,M*,p* kv ) 

Cn[*l,*2] ' ' ' 

respectively. 

Let ^1 = H a (l), < a ^ 1, and ^2 = L>i(B,(3) < (5 ^ 1. Let us construct an optimal with 
respect to order method for calculating integrals (7.1) on the classes of functions ^1 and surfaces ^2- 
Let S(G) be a "square" of the surface G. Divide the surface G into n parts g k , k = 1, 2, . . . , n, so that 
a "square" of each of the domains g k has the area of order \S(G)\/n, where | S'(G') | is the area of S(G). 
We take a point M k in each of domains g k at the center of the domain g k . 
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Calculate integral (7.1) by the formula 



Jf 



E/(M fe ) // 

k=l n . 



dS 



m ((ri - ti) 2 + (r 2 - t 2 y + (r 3 - t 3 y) 

Theorem 7.1. Formula (7.3), has the error 



+ Rn(f,G). 



(7.3) 



and is optimal with respect to order on the classes = H a , < a ^ 1, and ^ 2 = L±(B, (3) < /3 ^ 1 
among all formulas (7.2) with p = 0. 

Proof. Assume for simplicity that the surface G is given by the equation z = ip(x,y), (x,y) G 
Go, (f(x,y) ^ 0. Let (p x (x,y) := p, ip y (x,y) := q. Write the integral Jf as 



Jf 



f(n, T2, ^(ri,r 2 ))v / l +P 2 (ti,t 2 ) + q 2 {Ti,T 2 )dT 1 dT 2 



Go 



(ri - ti) 2 + (r 2 - t 2 ) 2 + Mri,r 2 ) - ip(h,t 2 )) 2 



A • 



(7.4) 



The function /(ti, r 2 , <^( T ii r 2 )) belongs to the Holder class #q, over Go, and the function 



TX 



is positive. 



(Tl-tl) 2 + (T2-t2) 2 +(v(-ri,T 2 )-^(tl,t2)) 2 

Let Mfc = (m^, m\,m\) be the nodes of cubature formula (7.2). Let i/)(t) := (d(r, {M^}))", where 
d(r,{Mk}) is the distance between the point r and the set of the nodes {M&}, where the distance 
is measured along the geodesies of the surface G. This distance satisfies the Holder condition H a (l). 
Hence the function 4>*{ti,t 2 ) = ^{t\, t 2 , tp(ri, r 2 )) belongs to the Holder class H a (A) and vanishes at 
the nodes (m\, m^), k = 1,2, ... ,n. Thus, 



Cn(*l,*2) > 



1 



S(G ) 



Go Go 



V>(ri,T2,y(ri,T 2 ))v / l +p 2 + q 2 dT l dr 2 dt 1 dt 2 
{n - h) 2 + (r 2 - t 2 ) 2 + (^(ri,r 2 ) - <^i, t 2 )) 2 



5(^o) / / ^ ri ' T2 '^ ri ' r2 )) c(Tl(ir2X 

Go 



x mm 

t 



^l+p 2 + q 2 



Go 



(ri - ti) 2 + (r 2 - t 2 ) 2 + Mn, r 2 ) - ¥>(*i, t 2 )) 2 

ff — 

J J (r(t,T 



-dTidT 2 ^ 



4 

^ — mm 



n; 



a/2 t 



■))' 



where S(Gq) is the "square" of the surface Go- 



Therefore the error of formula (7.3) is estimated by the inequality R n ^ r5~72- 
Theorem 7.1 is proved. ■ 

Remark 1. The method of decomposition of the domain G into smaller parts gu, k = 1,2, . . . ,n, 
described below, is optimal with respect to order for classes of functions ^1 = H a , < a ^ 1, and of 
surfaces * 2 = L (B, (3), < (3 < 1 for a < /?. 
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Remark 2. From formula (7.4) it follows that if the function f £ W r,r (l) and the surface G £ 
L s (B,a), then the function /(ti, T2, <p(ti, T2)) G W v ' v (A), where v = min(r, s). Therefore, repeating 
the above arguments, one proves that the accuracy of calculation of integral (7.4) by cubature formulas 
using n values of integrand function does not exceed 0(n~ v l 2 ). 

From this remark it follows that if the surface G consists of several parts, for example of surfaces 
G\ and G2 having common edge L, then it is necessary to calculate the integrals for the surface G\ 
and the surface G2 separately. If the surface G is divided into smaller parts gk, k = 1,2, . . . ,n, the 
domains gk, the curve L passes inside of these domains, should be associated with the class of surfaces 
Lq(B, 1). In these domains the accuracy of calculation of the integral does not exceed than 0(n^ 
where is the number of nodes of the cubature formula used in the domain g^ . 

For this reason the cusps and the nodes, in which three or more domains Gk, which are parts of the 
domain G touch each other, must belong to the boundaries of the covering domains gk, k = 1, 2, . . . , n. 

The universal code for computing the capacitances, described in Section 9, is based on optimal 
with respect to order cubature formulas for calculating integrals on the classes of functions H a , < 
a ^ 1 and of surfaces Lq(B,/3), B = const, a ^ (3, (3^1. 

The algorithm constructed in Section 9 is optimal on this class of surfaces and does not require 
special treatment of edges and conical points of the surface. 

When one studies cubature formulas on the classes W r ' r (A), r > 1, and L S (B, (3), s ^ 1, 0^/3^ 
1, one has to develop a method to compute accurately the integrals in a neighborhood of the above 
singular points of the surface. 

8. Calculation of weights of cubature formulas. 

In calculating weakly singular integrals by cubature formulas (6.8) it is necessary to calculate 
integrals of the form of 

dr\dT2 



Jkl{t\,t2) = j 



. ((n - hy + (r 2 - t 2 y) x 

for different values (ti,t2) S [—1, l] 2 . 

Let be (ii , ^2) £ \j- Let us consider two possibilities: 1) the square and the square Ajj have 
nonempty intersection; 2) the square A^; is does not have common points with the square Ay. 

First consider the second case, when the function 

((ti - ti) 1 + (r 2 - t 2 y) 

is smooth. Here (ti,t 2 ) G A^, and (ii,<2) € A^. 



In this case one has: 



r!2 2r 

< 



(( Tl _ tl )2 + (r2 _ t2) 2 ) A+r/2 
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and, if the squares A^ and Ajj do not have common points, one gets: 



9 r v?(ri,r 2 ) 



< 



2 r r\n 2X+r 



Similar estimates holds for partial derivative with respect to r 2 . 
Calculate the integral Jfcz(ii,i 2 ) by the Gauss cubature formula: 



((ri-ti) 2 + (r 2 -t 2 ) 



2\ x 



dTidT 2 + Rmm(Akl), 



where P m m = Pm^m > Pm (* = 1; 2) is the projection operator onto the set of interpolation polynomials 
of degree m with nodes at the zeros of the Legendre polynomial, which maps the segment [—1,1] onto 
the segment [a;fc,Xfc+i] for i = 1, and onto the segment [x/,x/+i] for i = 2. 



An integer m is chosen so that \R T 



< n 



-2-a 



for cubature formulas on the Holder class H n 



and |-R mm | < n r "for cubature formulas on the class W rr . 



This requirement is made because the error of calculation of the coefficients Jki(ti,t2) must not 
exceed the error of formula (6.5). 

Using r derivatives of the integrand in the error i? mm (A^), one gets: 



|-Rn»m(A W )| < 



B r 2 r r\ (2\ 
n J 



m 



r-l 



2-2A 



where B r is the constant appearing in Jackson's theorems. It is known that the constant B r are 
bounded by a constant, denoted b, uniformly with respect to r. In the case of periodic functions 6=1 
([18]), and in the general case b is apparently unknown. 

If r = 2 and m = B r 2 r r\n 2X , then one gets the error estimate given for cubature formula (6.5). 

Now, consider a method for calculating the integrals Jki{ti, t 2 ) when the square A^ has nonempty 
intersection with the square A^. For definiteness we consider the calculation of the integral Jjj(ti,i 2 ) 
by the formula: 



Jij(h,h) — J 



1 



((ri-ti) 2 + (r 2 -t 2 ) 2 ) A + / l 



dTidT 2 + R mm (Aij), 



where h = const > will be specified below. 
One has: 



\Rmm(\j)\ < h J 



dr\dT2 



'. (fa - ti) 2 + (r 2 - t 2 ) 2 ) X (((n - hY + (r 2 - t 2 ) 2 ) A + h) 



+ 



I 



A, 



(( Tl _ tl )2 + (r2 _ t2) 2 ) A + ^ 



dT\dT 2 = n + r 2 , 
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where D mm = I — P mm , and I is an identity operator, and 

<h [ dnd,T2 

Tl - ,,,,,,, / „ \(l+A)/2 - 

Ay ((ri - ti) 2 + (r 2 - t 2 ) 2 ) (A+1)/2 (((n - ti) 2 + (r 2 - t 2 ) 2 ) A + fc) 

The function — — — ^ — -r — is infinitely smooth. Using bounds for its first derivatives for 

A > 1/2, one gets: 

re 4 Arm 

From inequality (8.1) it follows that for getting accuracy 0{n~ l ~ a ) one has to have h = 
n -2(2A+a)/(i-A) an( j £ rom i ne q Ua iity (8.2) it follows that one has to have m = moi([nt 8A+4a ^ 1_A ' +a_3 ], 1). 

9. Iterative methods for calculating electrical capacitancies of conductors of arbitrary 
shapes. 

Numerical methods for solving electrostatic problems, in particular, calculating capacitancies of 
conductors of arbitrary shapes, are of practiacl interest in many applications. Electrostatic problems 
solvable in closed form are collected in [19,20,21]. Some of the problems were solved in closed form 
using integral equations, Wiener-Hopf and singular integral equations [22]. Electrostatistic problems 
for a finite circular hollow cylinder (tube) were studied in [23] by numerical methods. In [24] the 
variational methods of Ritz and Trefftz are discussed. Galerkin's and other projection methods are 
studied in [25]. In practice these methods are time-consuming and variational methods in three- 
dimensional static problems probably have some advantages over the grid method. There exists a 
vast literature on calculation of the capacitances of perfect conductors [20,26]. In [20] there is a 
reference section which gives the capacitance of the conductors of certain shapes. In [26] and [27] 
a systematic exposition of variational methods for estimation of the capacitances is given. In [28] 
there are some programs for calculating the two-dimensional static fields using integral equations 
method. In the monograph [7] iterative methods for solving interior and exterior boundary value 
problems in electrostatics are proposed and mathematically justified. Upper and lower estimates for 
some functionals of electrostatic fields are obtained. Such functionals are the capacitances of perfect 
conductors and the polarizability tensors of bodies of arbitrary shape. These bodies are described by 
their dielectric permittivity, magnetic permeability and conductivity. They can be homogeneous or 
flaky. The main point is: these bodies have arbitrary geometrical shapes. 

The methods, developed in [7], yield analytical formulas for calculation of the capacitances and 
polarizability tensors of bodies of arbitrary shapes with any given accuracy. Error estimates for 
these formulas are obtained in [7] . We give here the formulas for calculating the capacitances of the 
conductors of arbitrary shapes [7]: 
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\ -l 



(2k) 



dsdt 
r Jt r st Jr 



V>(Mi) . . . ip(t n -i,t n )dti ...dt n 



where S is the surface area of the surface T of the conductor, sq is the dielectric constant of the 

d 1 

dN t r a t ■ 



medium, r st := \s - t\, and ip(t,s) := gfr:r-, 



C (0) 

It is proved in [7] that 



AtteqS 2 



J 



<C, J 



dsdt 
rJr r st 



S = measr. 



c - c (n) 



< Aq n , 0<q<l, 



where A and q are constants which depend only on the geometry of F. 

We use these formulas are used to construct the computer code for calculating the capacitances 
of the conductors of arbitrary shapes. 

It is proved in [7], that 

-1 

r^6 n (t)dtd8 I , (9.1) 



= AtteoS 2 




where 8 n is defined by the iterative process: 



<Wi = —A6 n , S = 1, J 5 n dt = S, 
r 



(9.2) 



and A is defined by the formula: 



0N S 2irr st 
r 

where N s is the outer unit normal to V at the point s. 

To use iterative process (9.2), one has to calculate the weakly singular integral 

d 1 



dN s r st 



dt. 



(9.3) 



Let us describe the construction of a cubature formula for calculating integral (9.3), assuming for 
simplicity that the domain G, bounded by the surface T, is convex. This asumption can be removed. 

Let S be the inscribed in the conductor sphere of maximal radius r* , centered at the origin. 
Introduce the spherical coordinates system (r,4>,6), and the set of the nodes (r*, 9i), where <f>k = 
2/c7r/n, k = 0, 1, . . . , n, 61 = irl/m, I = 0, 1, . . . , m. Assume that m is even, and cover the sphere S 
with the spherical triangles k = 1, 2, . . . , N, N = 2n(m — 1). 

Let us describe the construction of the spherical triangles. For < < ir/m the triangles 
A fc , k = 1,2, . . . ,n have vertices (r*,0,0), (r* ,<f> k -i,0i), (r*,<f> k ,0i), k = 1,2, ... ,n. 
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For 9i < 9 < 0/+i, I = 1,2, ...,m/2 - 1, the triangles A fc ,fc = n + 2n(Z - 1) + j, 1 < 
j < 2n are constructed as follows. The rectangle [0, 27r; 9i, Oi+i] is covered with the squares 
Aki = [^fe, 4>k+i\ ®h 9i+i], k = 0, l,...,n — 1. Each of the squares Am is divived into two equal 
triangles Aj^ and A M , k = 0, 1, . . . , n — 1, I = 1,2,... , m/2 — 1. The spherical triangles A^ and 
A^, k = 0, 1, . . . , n — 1, Z = 1,2,..., m/2 — 1, are images of triangles A^ and A 2 kl on the sphere S 

As a result of these constructions the sphere § is covered with triangles A&, k = 1, 2, . . . , N. 

We draw the straight lines through the origin and vertices of the triangle A&, k = 1, 2, . . . , N. The 
points of intersection of these lines with the surface T are vertices of the triangle A&, k = 1, 2, . . . , N. 
As a result of these constructions the surface T is approximated by the surface Tn consisting of triangle 
Afc, /c = 1, 2, . . . , N, and integral (9.3) is approximated by the integral 

u{8) = l m wsh dt - (94) 

We fix each triangle fc = 1, 2, . . . , TV, and associate with it a point £ A^, fc = 1, 2, . . . , N, 
equidistant from the vertices of the triangle Afc, k = 1,2, . . . , N. We calculate integral (9.4) at the 
points Tfc, k = 1,2, ... , iV, by the cubature formulas constructed in paragraphs 5-7 for the Holder 
classes. After calculating the values U (r^), k = 1, 2, . . . , N by these cubature the integral 

C^ 1 ) = -4tt£ ^ fy y r-^U^dtds 
\n r N ) 

is calculated, where = Ufa) for t G A^, k = 1,2, ... ,7V, Sat is area of the surface T^, is 
approximation to the value of . The successive iterations are calculated analogously. 
10. Numerical examples. 

In this section the numerical results are given. As an example we calculated the capacitances of 
various ellipsoids, because for ellipsoids one knows ([19]) the analytical formula for the capacitance, 
which makes it possible to evaluate the accuracy of the numerical results. Consider the ellipsoid: 

x 2 y 2 z 2 
\- - — I = 1. 

a 2 b 2 c 2 

It is known [19,20] that the exact value of the capacitance of ellipsoid with a = b is: 



= 47rg 0V / (a 2 - c 2 ) 
arccos (c/a) 

Let a = b = 1, and £q = 1. We have calculate the capacitance C for different values of the 
semiaxis c. The results of the calculations are given in Table 1. 

It is known ([7], p. 43), that the capacitance of a metallic disc of radius a is C = 8aeo, and one can 
see from Table 1, that asymptotically, as c — > 0, this formula can be used practically for the ellipsoids 
with c < 0.001 with the error about 0.005. 



39 



Table 1 



c 


n 


m 


N 


Exact value 


Error 


Relative error 


Calculation time 


0.9 


40 


30 


2320 


12.144630 


-0.221200 


0.018212 


25 sec 


0.5 


40 


30 


2320 


10.392304 


-0.222042 


0.021366 


25 sec 


0.1 


40 


30 


2320 


8.5020638 


-0.301189 


0.035425 


25 sec 


0.01 


40 


30 


2320 


8.050854 


0.072132 


0.008959 


25 sec 


0.001 


40 


30 


2320 


8.005092 


-0.821528 


0.106374 


25 sec 


0.0001 


40 


30 


2320 


8.000509 


-1.068178 


0.133513 


25 sec 


0.9 


50 


40 


3900 


12.144630 


-0.180510 


0.014801 


1 min 15 sec 


0.5 


50 


40 


3900 


10.392304 


-0.185642 


0.017860 


1 min 15 sec 


0.1 


50 


40 


3900 


8.5020638 


-0.288628 


0.033947 


1 min 15 sec 


0.01 


50 


40 


3900 


8.050854 


-0.372047 


0.046212 


1 min 15 sec 


0.001 


50 


40 


3900 


8.005092 


-0.586733 


0.073295 


1 min 15 sec 


0.0001 


50 


40 


3900 


8.000509 


-0.933288 


0.116653 


1 min 15 sec 


0.9 


60 


50 


5880 


12.144630 


-0.152009 


0.012516 


4 min 


0.5 


60 


50 


5880 


10.392304 


-0.160023 


0.015391 


4 min 


0.1 


60 


50 


5880 


8.5020638 


-0.283364 


0.033328 


4 min 


0.01 


60 


50 


5880 


8.050854 


0.532250 


0.061110 


4 min 


0.001 


60 


50 


5880 


8.005092 


-0.391755 


0.048939 


4 min 


0.0001 


60 


50 


5880 


8.000509 


-0.880394 


0.110042 


4 min 
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